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ABSTRACT 

Recently Hall & Seddougui (1989) considered the secondary instability of large amplitude 
Gortler vortices in a growing boundary layer evolving into a three-dimensional flow with wavy 
vortex boundaries. They obtained a pair of coupled, linear ordinary differential equations 
for this instability which constituted an eigenproblem for the wavelength and frequency of 
this wavy mode. Investigations into the nonlinear version of this problem by Seddougui & 
Bassom have revealed several omissions in the numerical work of Hall & Seddougui; these 
issues are addressed in this note. In particular, we find that many neutrally stable modes are 
possible; we derive the properties of such modes in a high wavenumber limit and show that 
the combination of the results of Hall & Seddougui and our modifications lead to conclusions 
which are consistent with the available experimental observations. 


x Thi8 research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No.NASl-18605 while the authors were in residence at the Institute for Computer Applications in 
Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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§1. Introduction. 


The purpose of this article is to repeat and improve the calculation of Hall & Sed- 
dougui (1989) (hereafter referred to as HS), concerning an asymptotic description of 
the three-dimensional breakdown of steady, spanwise periodic Gortler vortices. These 
authors noted that in the experiments of Bippes & Gortler (1972) and of Aihara & 
Koyama (1981) this breakdown led to a time-periodic flow with wavy vortex bound- 
aries similar to those which occur in the Taylor problem. In order to investigate this 
phenomenon theoretically, HS superimposed small spanwise periodic travelling waves 
on the Gortler vortices and monitored their development. 

The first analytical work on Gortler vortices concentrated on the linear stability 
of external flows over concave walls. However Hall (1982a, b, 1983) showed that much 
of this early work was fundamentally flawed for it invoked the parallel flow approx- 
imation; but Hall demonstrated that this assumption is unjustifiable except in the 
limit of small vortex wavelength. Moreover, in this limit, the Gortler instability may 
be described by an asymptotic structure which accounts for boundary layer growth in 
a rational manner. This asymptotic structure was obtained by Hall (1982a) for the 
case of infinitesimal amplitude vortices and this was used by Hall (19826) to determine 
the modified structure in the case of weakly nonlinear vortices. Subsequently, Hall 
Lakin (1988) used this latter work to deduce the flow configuration for fully nonlinear, 
high wavenumber vortices, at which point the mean flow correction generated by the 
presence of the vortices becomes as large as the basic (undisturbed) flow itself. These 
fully nonlinear vortices are of the type whose stability to travelling wave disturbances 
were considered by HS. Hall & Lakin (1988) demonstrated that for these large ampli- 
tude vortices the flow structure consists of essentially three distinct regions. The main 
vortex activity is restricted to a central ‘core’ region which is bounded by two thin 
shear layers. The vortices decay exponentially within these shear regions and outside 
these zones the mean flow is governed by the usual boundary layer equations. 

Hall & Seddougui (1989) imposed infinitesimal secondary instabilities upon the 
flow within the shear layers; these modes took the form of short wavelength, high 
frequency travelling waves which were 7r/2 radians out of phase with the fundamental 
in the spanwise direction so that any instabilities which occurred produced locally 
wavy vortex boundaries in the shear layers. It was shown in HS that the governing 


1 


equations for these secondary modes take the forms 


d 2 u 
di 7 2 



2 i 2 2>/6 

rjv — —A v — V v H Kw = 0, 

3 3 


(l.la) 


d 2 w , _ v 9 i-v/6 

^^2 — ? (fit? + if) w + 2V w — — — — (fi77 + AT) Fv = 0. (1.16) 

In this pair of coupled ordinary differential equations for the functions 11(77) and w(rj), 
77 is an 0(1) co-ordinate based upon the thickness of the shear layers, (n,7n) are the 
normal and spanwise components of the velocity of the travelling wave disturbance, 
and AT and Cl are the (dimensionless) wavenumber and frequency of the imposed per- 
turbation. Further, the function V(tj) satisfies the Painleve equation, 


d 2 V 
d 7? 2 


-r]V = 



( 1 . 1 c) 


with V ~ as 77 — * —00 and V ~ \[2 At[rf) as 77 — > 00, see Hastings & McLeod 

(1980). For a complete description of the derivation of (1.1), together with an extended 
account of the previous theoretical and practical work relating to Gortler vortices, see 
HS. 

Hall & Seddougui (1989) concentrated on locating solutions of (1.1) for which the 
flow is neutrally stable, or, in other words, on finding solutions of (1.1) for which K 
and 12 are real. To ensure that the travelling Waves were confined within the shear 
layers, it was necessary to impose the boundary conditions 

v,w * 0 as 77 >±00. (l.lrf) 

This then implies that as 77 — ► 00, (l.la, 6) can be written as 



(1.2a) 


and 


d 2 w 

dr] 2 


i(Clr) + K)w = 0, 


(1.26) 


so that in this limit two independent solutions for v and w can be found in terms of the 
Airy function, At. For 77 — > —00, HS demonstrated that (l.la, 6) assume the forms 


d 2 v 2 iCl 2 iK 2\/6, i_ 

= < L3 “) 
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( 1 . 36 ) 


^ - (2 + in)rjw - iKw = + l -^I<{-ri)’v. 

dr) 1 3 o 

Then the appropriate expansions within (1.3) take the form 

v = (n O i+...)exp(-#7l0 , w = M’ (u>oi + • • •) exp (-<%!’) , 

where 

243^> 4 + 36</> 2 (6 + 5ift) - 32ft 2 = 0, 

and the two roots of this equation with positive real part were used in order to generate 

two independent solutions of (1.3) with v,w * 0, 77 > — oo. 

These asymptotic solutions for v and w as \r]\ — ► oo were taken as initial values 
in the numerical integration scheme used to solve (1.1). These equations were written 
as a system of four first order differential equations which was solved using a standard 
fourth order Runge-Kutta method. The integration procedure was started at 77 = —00 
and at 77 = 00 and continued to 77 = 0, finding two independent solutions from each 
direction. At 77 = 0 the continuity of a linear combination of the independent solutions 
from each direction produced a problem of the form 

Ax = 0, (1-4) 

where A (AT, ft) is a 4 x 4 complex valued matrix and x is a vector containing the 
coefficients of the independent solutions from 77 = ±00. Clearly there is only a non- 
trivial solution of (1.4) if det(A)= 0. HS found real values of I< and ft for which 
det(A)= 0 (and hence for which there are non-trivial neutrally stable travelling waves) 
by employing a Newton— Raphson iteration scheme for tw T o variables. 

Hall & Seddougui (1989) focussed their attention on the region AT, ft > 0 and 
found only one eigenvalue pair, which was located at 

{K, ft) = (4.156,0.742). (1-5) 

They speculated that other eigenpairs might possibly exist at higher values of It and 
ft, although none were found. Additionally, HS compared these theoretical findings 
with the results of the experimental observations by Kohama (1987) and Peerhossami 
Sz Wesfreid (1988). They reported good qualitative agreement with these practical 
results, although the lack of details given by these papers concerning the specific 
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experimental configurations used prevented a detailed quantitative comparison being 
made. 

Recent work by Seddougui & Bassom (1990) concerns with the extension of HS 
to the nonlinear regime. Seddougui & Bassom have shown that at the point at which 
the secondary wavy mode becomes nonlinear, the steady vortex flow of HS is affected 
by self-interactions of the wavy mode and the problem is then governed by an eighth 
order coupled set of ordinary differential equations. In the course of this work we 
reconsidered the linear problem of HS and discovered that the numerical results quoted 
in that paper are incomplete. In particular, we found eigenvalue pairs lower than (1.5), 
have shown that there are (plausibly) an infinite number of real valued solution pairs 
of (1.1), and have obtained an asymptotic description of the solution of this equation 
for K 1. 

The procedure for the remainder of this note is as follows. In section 2 we present 
a revised solution of (1.1) for 0(1) values of K and 12, in section 3 we consider the 
case K » 1 and finally we draw some conclusions. 


§2. The numerical solution of (1.1)— (1.3). 

To obtain real eigenvalues of the system (1.1)— (1.3) we employed the numerical 
method briefly described in the previous section with one modification. Instead of 
using the double Newton— Raphson iteration scheme to locate the neutral modes, we 
considered a whole range of real pairs (K, 12) and for each chosen pair we computed 
the matrix A defined by (1.4). Then, using a standard package, we constructed the 
contours in the ( K , 12) plane on which Re(detA)= 0 or Im(detA)= 0. These lines, for 
the region 0.05 < fi < 3, 0.1 < K < 7.9 are illustrated in Figure (1), where solid lines 
denote contours on which Re(detA)= 0 and broken lines those on which Im(detA)= 0. 
Plainly, where these contours cross we have a solution for infinitesimal, neutrally stable 
wavy modes. We remark that the apparent multi-crossings of the contours on the left 
hand side of this figure are not real features and are a direct result of the combination 
of very small absolute values of the determinant of A in this area and of the coarseness 
of the grid from which the contours were interpolated. By repeatedly refining this grid 
we can show that in reality there are no neutral pairs (K,£l) in this region. 

Having identified approximate values of possible wavenumbers and frequencies for 
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this secondary mode using this contouring technique, we then applied the Newton- 
Raphson iteration scheme to obtain more accurate values for these parameters. Over- 
all, we found eight eigenpairs, namely 

(0.690,0.372), (2.900,0.659), (4.156,0.742), (5.435,0.795), 

(7.53,1.00), (9.60,1.17), (11.4,1.27) and (15.7,1.60), 

of which the lowest five are marked on figure (1). We found no more eigenpairs 
within the size of grid used (0 < K < 16, 0 < ft < 3), but there is no reason to 
doubt that if this area were suitably enlarged more neutral pairs could be identified. 
Importantly, we deduce from (2.1) that there are neutral pairs lower than (1.5). In 
a practical setting, if the frequency of the imposed perturbations on a steady vortex 
flow were gradually increased from zero, then the mode with the lowest frequency 
eigenvalue would be expected to be the most dangerous as it would occur first. The 
eigenfunctions corresponding to this mode are shown in Figure (2) and the asymmetry 
of these functions is noticeable. Also, for the case plotted here, Re(u)> 0 across the 
whole of the dominant part of the mode, whereas the spanwise velocity component w 
has a much more oscillatory nature. 

For the higher eigenvalues it was found that it became increasingly more difficult 
to locate neutral modes accurately and although we have given only eight pairs in 
(2.1), we believe that there is indeed an infinite sequence of neutral modes. Inspection 
of the eigenfunctions of Figure (2); of those corresponding to the pair (4.156,0.742) 
(presented in HS); and of those corresponding to (9.60,1.17), illustrated in Figure (3), 
suggests some definite trends. In particular, the eigenfunctions are effectively confined 
to a region in space centred somewhere to the left of fj = 0. Furthermore, for the 
higher modes, the majority of the disturbance shifts to increasingly more negative 
values of r], and becomes ever more symmetrical in appearance. The ratio of the 
spanwise and normal velocity components, |w|/|v|, increases. These behaviours led 
us to an analytical consideration of the problem (1.1)-(1.3) in the limit K 1 in an 
attempt to elucidate the governing behaviour of neutral modes in this high wavenumber 
limit. This work is the subject of the following section. 


§3. The high wavenumber solution of (1.1)-(1.3). 
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Here we consider the solution of ( 1 . 1 )— (1.3) for K > 1 . It is found convenient to 
seek neutral modes in which the wavenumber K and the frequency ft take the forms 


k ={p + 0tKi + ...), n = n 0 P’, 


(3.1) 


where ^ > 1 and fto, Jli are 0(1) constants. Further, we suppose that the mode is 
concentrated within a region of thickness O(0~ * ) centred at a point 0(0*) from the 
origin r) = 0. If we then write 


K ■ i_ 

*1+ jj- = -P'ri c + 0 


(3.2a) 


and note that since for large negative rj, V ~ (— r/)a -f . . ., we also have 

rKy-fc) 


V = /3Nr,, + ±- 

r fto 


1 - 


2 + it ) 


+ 0(0 6 ) 


(3.26) 


These scalings suggest that 

v = VQ + 0~ r Vl + . . . , W = 0* (wQ + 0~ *Wl + . . . ) , 


(3.3a) 


and inserting these expansions in ( 1 . 1 ) and equating coefficients of leading powers of 0 
we find that r\ c = 0. At next order we obtain the coupled equations for the functions 
vq and wq, 


d 2 vo 2 ifto 2\/6 

yv o + - t== wq = o, 


dy 2 3 ' 3Vft^ 

d 2 wo . ( 2 o 


dy 2 

with boundary conditions 


+ ( - mow) m-- 


yv 0 = o, 


(3.4a) 

(3.46) 


^ 0,^0 >0 as y ► ±oo. ( 3 . 4 c) 

By eliminating wq between (3.4a, 6 ) we can derive a fourth order differential equa- 
tion for no which may be solved numerically using standard Runge-Kutta ideas. How- 
ever, careful inspection and analysis of (3.4) reveals that a solution of the system for 
which vo, wq — > 0 as J 77 I — > 00 is possible for all non-zero real fto. Without loss 
of generality, we can find solutions of (3.4) such that Re(vo) and Re(zco) are even 
functions about y = 0 whilst Im(vo) and Im(tvo) are odd functions. As fto is not 
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determined at this stage, we are forced to consider the next order equations. These 
are found to be 


d 2 v\ 2?fio 2\/6 V6fio f -KjA 

W - — yn + wrr 0 m = — y y -T ^) W0 


(3.5a) 


d z w i 

dy 2 




W 1 


iV^O 


-yv i = 



(3.56) 


with 


vi,wi — ► 0 as y — > ±oo. 


(3.5c) 


Now as the homogeneous forms of (3.5) are precisely those equations given in 
(3.4), the set (3.5) has a solution only if a certain compatibility condition holds. To 
derive this compatibility criterion we consider the system adjoint to (3.4) which in this 
case is formed by the functions (yi(y), 92{y)) where 


with 


d^g\ 2iflo i\/6fio 

vn-— r -vn = *. 

w + ^M 91 + \fr 0 - ,n ° y ) g2=0 ’ 


91, 92 ► 0 as y > ±oo. 


(3.6a) 

(3.66) 

(3.6c) 


As before, we can obtain solutions of (3.6) in which the real parts g\ and Q 2 are 
even functions whilst the respective imaginary parts are odd functions. We multiply 
(3.5a) by <71 , (3.56) by g 2 , add the resulting equations and integrate by parts over the 
range (—00, 00). Then on comparing imaginary parts of the resulting eigenrelation we 
find that (3.5) has a solution only if 


\/6f2o 


3 

/ oo r 00 yo° 

(j/Im(<7iico)) dy + 2 J ( yIm(g 2 WQ )) dy - -jL J ^ ( y 2 Re(g 2 Vo )) dy = 0, 


(3.7) 

where the evenness and oddness of the various functions about y = 0 has been utilised. 
Incidentally, we note that on equating real parts of the eigenrelation we obtain K\ = 
0. Numerically, we found that (3.7) is satisfied for flo ~ 0.26 and the respective 
disturbance functions do and w 0 for this case are presented in Figure (4). 

This then suggests that high frequency, high wavenumber modes are indeed pos- 
sible; in particular if K 1 we have ~ 0.2 6K » . Additionally, the position of the 
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disturbance within the shear layer moves towards the core region of vortex activity, 
the perturbation is confined within a thin, 0(/i“ »■) sized zone relative to the depth 
of the shear layer, and the spanwise velocity component of the wavy mode is 0(K‘) 
times that of the normal component. These asymptotic results provide satisfactory 
agreement with the higher numerical values calculated in (2.1) and it is observed that 
there is resemblance between the asymptotic eigenfunctions of Figure (4) and those 
corresponding to the eigenpair (9.60,1.17) shown in Figure (3). We also remark that 
here we have not specified K and have only assumed that this is some large quantity. 
In practice the actual possible values for K 1 would be determined by a higher 
order problem but we do not pursue such a consideration here. Instead, we feel that 
the principal result of this analysis is that, almost certainly, there is an infinity of 
neutral wavy modes for this flow structure. This confirms the conjecture presented 
by HS, who postulated that neutral wavy modes may possibly exist at wavenumbers 
higher than that for the mode with eigenvalues given by (1.5). 


§4. Conclusions. 

In this article, which should be read in conjunction with HS, we have demonstrated 
that there are many solutions for infinitesimal, neutrally stable wavy modes confined 
to the shear layers of a fully nonlinear, high wavenumber Gortler vortex. This is 
at variance with the result of Hall Sz Seddougui (1989) who found only one such 
mode: although their conclusion, that such a large amplitude vortex is unstable to 
modes trapped within the shear layers, remains unaltered. Our results are improved 
from those of HS due to the use of the contouring method described in section 2. 
This ensured that all possible neutral mode solutions within the region 0 < K < 16, 
0 < ft < 3 should be identified. This is in contrast to the scheme employed by HS, who 
used the Newton-Raphson technique alone and relied on finding all possible neutral 
solutions by using a sufficiently wide variety of initial estimates to which the iteration 
procedure was applied. 

Since the wavy modes are stable when their non-dimensional frequency Q, — 0, 
HS noted that as ft is increased the mode with the lowest frequency is potentially 
more dange ro us than the other modes as it occurs first. How'ever, at large values of 
ft it is unclear as to which is the most important mode as we would need to identify 
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the mode with the largest growth rate. The contour plot in Figure (1) shows that 
(K,n) = (0.690,0.372) is the lowest mode within the region IC > 0,0 > 0 and that, 
plausibly, an indefinite number of neutral modes exist. In particular there are modes 
with K » O > 1 whose asymptotic structure has been discussed in section 3. 

The physical and experimental implications of our revised analysis follow similar 
lines to those given in HS. Rather than repeat the discussion given there and to which 
the reader is referred, we merely emphasise that our results concerning the possibility of 
wavy modes in Gortler vortices are consistent with the available experimental accounts 
on the subject, although it is very difficult to obtain quantitive comparisons between 
our theory and the practical observations. 

Our conclusions concerning the existence of many possible neutral linear modes 
raises questions of which of these modes is likely to be the most important in practice. 
This may be resolved by pursuing weakly nonlinear and, ultimately, fully nonlinear 
analyses of these wavy modes — a task currently being tackled by Seddougui & Bas- 
som. When concrete conclusions become available from this study, the full implications 
of our work for this wavy vortex problem should become clearer. 
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FIGURE CAPTIONS 


Figure (1). Sketches of the contours in (K,Sl) space on which the real and imagi- 
nary parts of det(A), defined by (1.4), vanish. On solid lines Re(det A)=0, on broken 
lines Im(det A)=0. We have a solution for neutral wavy modes "wherever these con- 
tours meet and five such locations are ringed, see (2.1). Here 0.05 < ft < 3, and (a) 
0.1 < K < 4, ( b ) 4 < K < 7.9. 

Figure (2). Neutral eigenfunctions 17 ( 77 ) and w{q) for (K, ft) = (0.690,0.372). a) 
Re(v), b) Im(u), c) Re(u;), d) Im(u>). 

Figure (3). Neutral eigenfunctions 17 ( 77 ) and w(tj) for (K,fl) = (9.60,1.17). a) 
Re(v), b) Im(u), c) Re(ie), d) Im(w). 

Figure (4). Asymptotic eigenfunctions v 0 (y) and w 0 (y) (defined by (3.3)) in the 
limit K > 1 with flo = 0.26. 
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Figure 4(b) 
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